DISTRIBUTIONAL LIMITS FOR THE 
SYMMETRIC EXCLUSION PROCESS 



Thomas M. Liggett 
University of California, Los Angeles 

October 18, 2007 

Abstract. Strong negative dependence properties have recently been proved for the symmetric 
exclusion process. In this paper, we apply these results to prove convergence to the Poisson and 
normal distributions for various functionals of the process. 



1. Introduction 

The symmetric exclusion process on the countable set 5' is the Markov process rjt on {0, 1}"^ 
with formal generator 

r]{x) = l,r]{y)=0 

where rjo^^y is the configuration obtained from r] by interchanging the coordinates r7(x) and 
77(1/). Here p{x,y) = p{y,x) are the transition probabilities for a symmetric, irreducible, 
Markov chain on S. For background on this process, see Chapter VIII of Liggett (1985). 

Many limit theorems of various types have been proved for this process. Examples are the 
central limit theorems for a tagged particle and for the flux in one dimensional systems by 
Arratia (1983), Kipnis and Varadhan (1986), De Masi and Ferrari (2002), Jara and Landim 
(2006), and Peligrad and Sethuraman (2008). In this paper, we focus on limit theorems that 
can now be proved using the recently obtained strong negative dependence properties of the 
symmetric exclusion process. 

A probability measure n on {0, 1}"^ is said to be negatively associated if 

fgd^ < / fd^ / gdfx 
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for all increasing continuous functions /, g on {0, 1}*^ that depend on disjoint sets of co- 
ordinates. Theorem 5.2 of Borcea, Branden and Liggett (2008) asserts that if the initial 
distribution of the symmetric exclusion process r]t is a product measure, then the distribution 
of r]t is negatively associated for all t > 0. In fact, by Proposition 5.1 of that paper, it has a 
stronger and even more useful property, known as strong Rayleigh. 

Limit theorems for negatively associated random variables have been proved by a number 
of authors - see Barbour, Hoist and Janson (1992), Newman (1984) and Roussas (1994), 
for example. In the case of convergence to the normal law, none of these results quite fit 
our setting. In our situation, there is generally no translation invariance in the covariance 
structure, and the sum of off-diagonal covariances is often not "little o" of the sum of variances. 
However, we will see in Section 2 that the strong Rayleigh property makes it quite easy to prove 
convergence to the Poisson and Gaussian laws, given estimates of variances and covariances. 
Therefore, we will not need to use these earlier results. 

The first situation we will consider involves the extremal stationary distributions for the 
process. We recall their description - see Chapter VIII of Liggett (1985). Let 



and for a E Ti, let be the product measure with marginals u{r] : r]{x) = 1} = a{x). Then 
the limiting distribution ast ^ oo of the process at time t exists if the initial distribution is Ua, 
call it Ha- The result is that the extremal stationary distributions are exactly {//q, a e H}. If 
a is constant, then //q, = ly^ so we are really interested in nonconstant a's, in which case very 
little is known about the corresponding stationary distributions other than the marginals - 
Haiv '■ r]{x) = 1} = a{x). If p{x, y) are the transition probabilities for simple random walk on 
a homogeneous tree, for example, there are many such nonconstant ct's, and therefore there 
are many extremal stationary distributions that are not of product form. We now know from 
the results in Borcea, Branden and Liggett (2008) that is negatively associated - and in 
fact strong Rayleigh - for each a E Ti. 

We will use the following notation. For n > 1, p^'^\x,y) are the n-step transition proba- 
bilities for the Markov chain with transition probabilities p{x, y), and for t > 0, 



are the transition probabilities for the corresponding continuous time chain Xf. The Green 
function is given by 




y 





The Dirichlet sum of an o: G is defined by 



x,y 
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This quantity is finite for many, but not all, elements of 7i if S' is a regular tree, for example. 
A construction of a class of infinite graphs with only one end that support nonconstant a & H 
with ^{a) < cxD is constructed in Cartwright and Woess (1992). In this context, we have the 
following results. We use =^ to denote convergence in distribution. 

Theorem 1. Suppose a ETi and $(a) < oo. If <Z S satisfy 
(a) limn_^oo supa;£5„ a{x) = 0, lim„_^oo Ea;eS„ ^i^) = A < oo, 
and 

(h) sup^^„ Ej/6S„ G{x, y) < oo, 
then under /la, 

r]{x) ^ Poisson{\). 

Theorem 2. Suppose a E H and $(a) < oo. If Sn C S satisfy 

lim Q:(a;)[l — a{x)] = oo 

n—^oo ' 

and 



(1.1) liminf ^-%^ > $(a) 

n-oo sup^ Eyes„ Gix, y) 



then under Ha, 



Furthermore, 



Ea:eS„ Vi^) T^xeSr, Q'(^) j^(r. -jN 



(1.2) Far^„ J] 77(0:)) ^ Yl " «(^)]' 

(which is at most one) is bounded below by a positive constant. If the left side of (1-1) is 
infinite, then the limit of (1.2) as n — > 00 is 1. 

Theorems 1 and 2 will be proved in Section 3, after deriving limit theorems for general 
strong Rayleigh Bernoulli random variables in Section 2. 

As an example of the application of Theorems 1 and 2, let S be the binary tree, and 
let the chain have nearest neighbor jumps with probability 1/3 each. Write S = L U R 
where L, R are defined as follows: An basis edge is fixed, and its endpoints are called left 
and right respectively. Then L is the set of vertices that are closer to the left vertex than 
to the right, and R = S\L. Each vertex a; in S' is assigned a level l{x) > 0, which is 
the distance from x to the closer of the two endpoints of the basis edge. Then G{x, y) = 
2-'^(^.y)+i^ where d{x, y) is the distance between x and y. Therefore, sup^ 'Yliyi{y)<n ^i^i v) ~ 

3n, SUp^ J2y&L:l{y)<n V) = '^'^ ^ud SUp^ J2yeL:l{y)=n y) = 3 - 2" 



-n 



4 



THOMAS M. LIGGETT 



For < A, p < 1, let a en be defined by 

a{x) — 

and put iJ, = fia- Then ^{a) = 2{p - A)^/9, 



^ + 3^ iixeL, 
P + 3W if a; e 



E'^ ^ ?7(x) = (A + p)(2"-l), and 

x:l{x)<n 

J2 Var^vix) = (2- - l)[p(l - p) + A(l - A)] + (A - pf 

x:l{x)<n 



2 4 

3" -5(1 -2-") 



It follows that for Sn = {x : l{x) < n}, the left side of (1.1) is infinite if 

p(l-p) + A(l- A) > 0. 
Therefore, for all choices of A, p, Theorem 2 implies that 

^x-l(x)<nV{x) - (A + p)2'^ 

(1.3) ^x.i{x)<n lyj ^ ^ jv(0,p(l-p) + A(l-A)). 

Next take A = 0, p = 1, in which case (1.3) has little content. Then Theorem 1 gives 

r]{x) =^ Poisson{l/3). 

x€.L:l{x)=n 

li Sn = {x e L : l{x) < n}, the left and right sides of (1.1) are | and | respectively, so (1.1) 
does not hold in this case. Nevertheless, we will see at the end of Section 3 that more careful 
estimates imply that 

(1.4) E.gL:K.)<n^(^)-i ^^^Q^^^^ 

where cr^/n is asymptotically between ^ and |. 

The next situation we consider was proposed by Pemantle (2000) as an application of 
the then hoped for negative dependence properties of the symmetric exclusion process. Now 
S = Z^, and p{x, y) — p{y — x), with |a;|p(a;) < 00. For the initial configuration, we take 

r 1 ifx<o, 

= { ifx>0. 

Since ^^^o^yP{x,y) <cx),Wt = Y.x>o^ti^) < °° ^- 
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Theorem 3. Suppose = Y^^n'^pin) < oo. Then 

Wt - EWt 
[Far(Wt)]V2 



Wt - EWt 

(1-5) r ^mi) 



as t ^ oo. Furthermore, 



and 



EWt o- 

lim 



tV2 

is bounded above and below by positive constants. 

Theorem 3 will be proved in Section 4. It seems likely that if the distribution p{-) is in the 
domain of normal attraction of a (symmetric) stable law of index a e (1, 2), then (1.5) holds 
with Var{Wt) of order t^/", but this has not been checked. 

The major result from Borcea, Brahden and Liggett (2008) that we use in this paper is 
Proposition 5.1, which asserts that the strong Raylcigh property is preserved by the evolution 
of a symmetric exclusion process. The proof given there uses results from earlier papers, 
including Obrcschkoff 's Theorem, which are not part of the toolkit of a typical probabilist. In 
Section 5, we present an elementary proof of that result in order to make the present paper 
essentially self-contained. 

2. Limit theorems for strong Rayleigh measures 
Consider a probability measure /x on {0, 1}"^. Its generating polynomial is defined for z & 



Mi) 



I 



by 

n 

Q{z) = Q{z^,...,Zn)^E^\[ 
The measure /x is said to be Rayleigh if 

for all z G -R", and to be strong Rayleigh if (2.1) holds for all z G R^. Note that when 
z = (1, 1, 1), (2.1) says that r]{i) and rj{j) are negatively correlated under n. By Theorem 
4.9 of Branden, Borcea and Liggett (2008), the strong Rayleigh property implies negative 
association. 

Branden (2007) proved that the strong Rayleigh property is equivalent to the following 
property, which is known as stability: Q{z) ^ if Zi has strictly positive imaginary part for 
each i. This is the key to the following representation. 
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Proposition 2.1. Suppose /j, is strong Rayleigh. Then there exist independent Bernoulli 
random variables d with parameters pi so that the distribution ofY^^ r}{i) under jj, is the same 
as that of^-Ci. 

Proof. Setting Zi = w in the expression for Q, we see that the polynomial in one variable 

Q*{w) = Q{w, w, w) = Ef'w^i "^'^ 

has no roots with positive imaginary part, and therefore all of its roots are real. Since 
Q*{1) — 1, it follows that Q* can be written in the form 

Q*{w) = l[{w + Wi) j ^(1 + Wi) 

i ' i 

where the Wi are real. Since Q*{—Wi) = 0, we see that Wi > for each i. Letting pi = 
1/(1 + Wi), this becomes 

Q*{w)^l[[piW + {l-pi)], 

i 

which is the generating polynomial for Yli d- 

Using this result, it is easy to extend the classical limit theorems to triangular arrays of 
strong Rayleigh random variables. 

Proposition 2.2. Suppose the Bernoulli random variables {rin{x)} are strong Rayleigh for 
each n. 

(a) //limn^oo Y.X ^Vnix) = A, limn^oo Y.xi^Vnix)]'^ = 0, and 

lim y^Cov{r]n{x),r]n{y))^0, 

xj^y 

then 

'^^Vnix) ^ Poisson{X). 



(b) Iflim.n^ooVar{Y,xVn{x)) = oo, then 



VVar{^^ r]n{x)) 



iV(0,l). 



Proof. Using Proposition 2.1, let be Bernoulli random variables that are independent in 
i for each n, and have the property that 



^r)n{x) and J^Cn,i 

X i 
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have the same distribution for each n. It suffices to show that the conditions for Poisson and 
normal convergence hold for the array Cn,i- But this is immediate from the assumptions and 
the following identities: 

X i 

and 

Vari ^r]n{x) \ ^^Var{r]n{x)) + ^Cov{r]n{x),r]n{y)) 

^ X ^ X x^y 

^ i ' i 

Remark. Proposition 2.2(a) is a consequence of Theorem 11 in Newman (1984) or Theorem 
3D in Barbour, Hoist and Janson (1992). The latter also gives bounds on the total variation 
distance from the Poisson distribution. We note that under the strong Rayleigh assumption 
we are making, the proof is very simple. 

3. The stationary distributions 

We begin the proofs of Theorems 1 and 2 by obtaining a bound on the covariances for 
the measure ji = jia- Let U and U{t) be the generator and semigroup for the motion of two 
independent copies of the Markov chain with transition probabilities pt{x, y), and V and V{t) 
be the generator and semigroup for the motion of two copies of that Markov chain with the 
exclusion interaction. Then by (1-28) and (1-29) on page 373 of Liggett (1985), 

Ef'r]{x)r]{y) = lim V{t)f{x, y), x y, 

t — >oo 

where f{x,y) = a{x)a{y). Since U{t)f{x,y) = f{x,y), the integration by parts formula gives 

-Cov^{rj{x),rj{y)) =f{x,y) - lim V{t)f{x,y) 

t— >oo 

= hm \u{t)f{x,y)-V{t)f{x,y) 

t—^00 

(3.1) ft 

= lim / V{s)[U -V]U{t-s)f{x,y)ds 

) 

V{s)[U-V]f{x,y)ds 

for X y. From page 366 of Liggett (1985), we see that for x ^ y, 

(3.2) {U - V)f{x, y) = p{x, y)[f{x, x) + f{y, y) - 2f{x, y)] = p{x, y)[a{x) - a{y)f. 



Define A{x,y) = p{x,y)[a{x) — a{y)]'^ for all x,y. 



8 



THOMAS M. LIGGETT 



Let gn{x,y) — Is^xSni^^v)- This is a positive definite function. Therefore, using the 
symmetry of V{s)^ the fact that A(a;, y) = if a; = y, and Proposition 1.7 on page 366 of 
Liggett (1985), we see that 

^ CoVi^{r]{x),r]{y)) = ^ gn{x,y)V{s)A{x,y)ds 

j-OO 

I ^^i^^yWi^)9n{x,y)ds 
^0 x,y 

^ A{x, y)U{s)gn{x, y)ds 



(3.3) "^'"^ 



x,y 



POO 

J2 A(x, y) / P^X, G Sn)PHX, G Sn)ds, 

x,y 



where Xg is the Markov chain with transition probabihties Ps{x, y). 

Proof of Theorem 1. Given the strong Rayleigh property of and Proposition 2.2(a), it 
suffices to check that 

hm V Cov^{r]{x),r]{y)) ^Q, 

x,y£Sn;xi^y 

and therefore we need to check that the right side of (3.3) tends to zero. To do so, note that 

/"OO 

(3.4) / P^(X, G Sn)ds = G{x,u) < Ci 

for some constant Ci by assumption (b). Next, since 

supe 77 = e < -= 
for some constant C2, we have for any N >1, 

py{x,eSr.)<j2J2p'''M + ^^ 

SO that (using the first part of assumption (a)) for each y 

(3.5) hm sup py(Xs G Sn) = 0. 

n-i-oo s>o 

It foUows that the right side of (3.3) tends to zero by (3.4), (3.5), the finiteness of $(«) and 
the dominated convergence theorem. 
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Proof of Theorem 2. By Proposition 2.2(b), we need to show that 
(3.6) hm Varai 'r]{x) ] = oo. 

By (3.3), 



CoVi^{r]{x),r]{y)) <^{a)sup ^ G{x,u). 

Therefore, (3.6) follows from the assumptions in the theorem. 

We now return to the example discussed in the introduction - the simple random walk on 
the binary tree with 

if a; e L, 
-'- ~ 3.2i(^) ii X & R, 



a{x) — 



and Sn — {x E L : l{x) < n}. We will see that Proposition 2.2(b) applies, even though (1.1) 
fails. In order to do so, we will use the structure of the problem to estimate the right side of 
(3.3) more carefully. 

First note that l{Xt) is a Markov chain on the nonnegative integers with drift i. Therefore 

(3.7) -r^3 

by the strong law of large numbers. Secondly, if 7 > are chosen so that 

e-^ + e+^ = 3(l-/?), 

then 

g/3t-7/(^t) 

is a supermartingale, so that 

1 > £;^e^''*-^^(^"*) > e^"*-^"P^(Z(X„i) < n). 

This implies 

P^'iliXnt) <n)< e""^^-^'^ < e^-^* 
if n > 1 and t > 7//?, which provides the domination in the following computation. 

1 r°° 

lim - V A(a;,y) / P^(X, e >Sn)P^(X, e Sn)ds 

x,y 

/•oo 

(3.8) = lim TA{x,y) / P^(X,i e Sn)PnXnt e Sn)dt 

x,y 

=3j]A(a;,y)[l-«(a;)][l-a(y)]. 



x,y 
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For the final step above, note that the integrand in the middle line tends to zero if t > 3 by 
(3.7), while if t < 3, 

lim eS^)=P''{X,eL eventually) = \-a{x). 

n— >oo 

Next, we compute the right side of (3.8). 

3^A(x,y)[l-a(a;)][l-a(y)]= J] [a{x) - a{y)f[l - a{x)][l - a{y)] 

x,y d{x,y)^l 

34 (-3 . 2n)2 3 . 2^ 3 • 2"-i ^ (3 • 2^)2 V 3 • 2^ A 3 • 2^-^ / 

40 1 
< -. 



189 3 

Combining this with (3.3) and (3.8), we see that 

liminf — > 0, 

n—K30 fi 

SO that Proposition 2.2(b) gives (1.4). 

4. Limit theorems in one dimension 

In this section, we will prove Theorem 3. We need to consider the first and second moments 
of Wt- By duality, 

EWt =E^J2^t{x) = Y,P^{'nt{x) = 1) 

x>0 x>0 

=^p-(x,<o)= Yl p^^'^^y) 

x>0 y<0<x 



(4.1) 



J2npt{0,n) = E^X+. 



n=l 



Similarly, 



(4.2) J][P^(r/,(a;) = l)]'= ^ Mx,u)Mx,v) = E^^''^ mmiX+ ,Y+), 

x>0 x>0;u,v<0 

where Yt is an independent copy of Xf. 

For the covariances, we proceed as in Section 3, but this time with f{x,y) = i{x,y<o} and 
9{x,y) = l{^,y>o}- For X y, {U - V)f{x,y) = p{x,y)[f{x,x) + f{y,y) - 2f{x,y)], so for 
xj^y, 

{U - V)U{s)f{x, y) =p(x, y) [U{s)f{x, x) + U{s)f{y, y) - 2U{s)f{x, y)] 
=p{x,y)[P^{X,<0)-py{X,<0)]^ 
=p{x,y)[P%X, >x)- P^X, > y)]2 = A,{x,y), 
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where we define for all x,y, 

p{x,y)[P^{x<Xs<y)f ifa;<y, 

As(x,2/) = <( p{x,y)[P'^{y<Xs<x)f if y < a;, 

ii y = X. 



Then letting 



we see that 



K{t) = - CovMx)My)), 

x,y>0;xj^y 



(4.3) 



K{t)^ J2 [U{t)-V{t)]f{x,y)^ Yl fv{t-s){U-V)U{s)f{x,y)ds 

x,y>0;x^y x,y>0;xj^y ° 

= / y23(^^y)^i^- ^)^s{x,y)ds= y2^s{x,y)V{t- s)g{x,y)ds 

x^y x,y 

■^0 x,y 

= I y2^s{x,y)P\Xt-s < x)P\Xt-s < y)ds. 
Jo ^ „ 



x,y 

Let p{t, x) = P^{Xt < x). Since As{x, y) = As{y, x) = As{l — x, 1 — y) , we may continue 
(4.3) by writing 



1 1 x 

K{t) <- V As(x, y)[p{t - s, x)p{t -s,y)+ p{t - s, 1 - x)p{t - s, 1 - y)]ds 
^•^0 x,y 

= / J]p(a;,y)[P°(a; <2/)]27(t-s,x,2/)ds, 



(4.4) 



x<y 

where 

lit, X, y) = p{t, x)p{t, y) + p{t, 1 - x)p{t, l-y). 
Note that since p(t, a;) + p(t, 1 — x) = 1, it follows that 7(t, y) < 1. Now let 

(4.5) r(t, n, tt, v)= Y^ '-^{t^x^x -\- n) < {n — \v — u\)'^ . 

x:x<u,v<x+n 

Then using the symmetry and translation invariance of Ps(a;, y), 

K{t)< / ^p(n)^ps(0,'u)ps(0,'i;)r(t- 5,71,^,^)^5 

n=l u,v 



(4.6) 



= / Yp{n)E''V{t-s,n,Xs,Xs-X2s)ds 

•^0 n=l 

< / Yp{n)E\n-\X2s\)^ds. 

Jo n=l 
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Now assume that cr^ = "'^^'("') ^ write 

n 

- \X2s\)+ = 5^ - |A;|)p2.(0, k) < n^p2s{0, 0), 

fc=— n 

where the inequality comes from 

r 1 ^ 

\P2s{o,k)f = 5;]p.(o,j)p,o-,A;) < J]b.(o,i)]'E[p-(^'^)]' = b2.(0,0)]^. 

Therefore, 

(4.7) hmsup ^ / y p(n)£;°(n - |X2s|)+ds < V n2p(n). 



n>N 



by the local central limit theorem. For the terms corresponding to small n, we have given up 
too much in using the inequality in (4.5). To handle these terms, first note that for fixed k, 
conditionally on \X2s\ = k, Xs/^/s is asymptotically normal with mean and variance 
Now take fixed k, n with < k < n. Then 



£;0[r(t - s, n, Xs, Xs - X2s) I X2s = k] =£;° 



j{t — s,x,x + n) 



Xs —n<x<Xa —k 
n-1 



X2s = k 



^ 7(t - s,Xs-y,Xs-y + n) 

y=k 



X2s = k 



Write 



£;0[7(t -s,Xs- y,Xs -y + n)\X2s = k]= P'^iYts < X^ - y, Zts < X, - y + n \ X2s = k) 
+ P'^iYt-s <l-Xs + y,Zt-s<l-Xs + y-n\X2s = k), 

where X^, 1^, are independendent copies of the random walk. If s,t ^ oo with s/r ^ r & 
(0, 1), then the above expression converges to 



P N2 < 



r / r 

NuNs<J— -Ni 



+P{ N2 < 



2(1 -r) 



2(1 - r 



2(1 - r)' 



-NuN3<- 



2(1 - r 



where A^i, A^2, -^3 are independent normally distributed random variables with mean zero and 
variance 1. Call this expression h{r). Note that h{r) < 1 for r < 1, so that 



Jo 
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Passing to the limit in (4.6), we see that 



K{t) ^ 2 . ^ 1 V- 2 . X 

limsup — r=- < — r= > n pin) -\ r= > n pin). 

^ * n=l * n>iV 



Taking N ^ oo gives 



Iim sup — ^ < 



RecaUing from (4.1) and (4.2) that 

Var{Wt) = E^X+ - min(X+, Y+) - K{t), 

we see that 

hmsupt"^/Var(W^t) < 
t^oo 2y7r 

and 

liminf t-^/^Var{Wt) > '^^"^^^ > 
t—Kx> 2y7r 

So, Theorem 3 follows from Proposition 2.2(b) and the strong Rayleigh property of {r]t{x), x > 
0}. 

5. Stability and the Symmetric Exclusion Process 

In this section, we present an elementary proof of the basic fact needed to prove preservation 
of stability by the symmetric exclusion process. A similar proof was obtained independently 
by Borcea and Branden (2008), 

A polynomial Q{zi, Zn) with complex coefficients is said to be stable if Q{zi, Zn) ^ 
whenever 9(zi) > for each i. The key result needed to show that the generating polynomial 
of the distribution a symmetric exclusion process at time t is stable whenever this is the case 
at time is the following. 

Theorem 5.1. Suppose the multi-affine polynomial Q is stable. Then so is the polynomial 
Qp for < p < 1, where 

Qp{zi, Zn) = pQ{zi, Zn) + (1 " p)Q{z2, Zi, Z3, Zn). 



The proof of Theorem 5.1 is based on the following characterization of stability for multi- 
affine polynomials in two variables. 

Proposition 5.2. Suppose h{z, w) — a + bz + cw + dzw, where a, 6, c, d are complex, and not 
all zero. Then h is stable if and only if 

(5.1) 3?(6c - ad) > \bc - ad\, Q{ab) > 0, Q{ac) > 0, %bd) > 0, ^{cd) > 0. 
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Proof, li b = c = d = 0, the h is automatically stable, since then a ^ 0. If d = 0, 6 0, then 
h{z,w) = OiS 

(a + cw)b 



z = 



|2 



so stability is equivalent to Q{w) > ^ $5(2;) < 0, where 

•^(^) - j^j2 • 

Therefore, since ^{w) is arbitrary, stability is equivalent to c6 > and Q{ab) > 0, and these 
imply Q'(ac) = bcQ{ab)/\b\'^ > in this case. 

So, we may now assume that d ^ 0. Solving for z, we see that h{z, = iff 

(5.2) b + dw = 0, a + cw = 
or 

, , , a + cw ab + cbw + adw + cd\w\'^ 

(5.3) b + dwf^O, z— — 



(5.4) '"^ ~ ~ J ~ ~TjT2' bc=ad. 



b + dw \b + dw\'^ 

Case (5.2) occurs iff 

b _ _bd 

d~ \d\' 

Therefore, if be ^ ad, stability is equivalent to the statement 

(5.5) %w) > ^ '^{ab) + ^{w)'^{cb + ad) + '^{w)^{cb - ad) + \w\'^'^{cd) > 0, 

while if be = ad, stability is equivalent to this, together with '^{bd) > 0. 
Letting w be real, we see that for (5.5) to hold, we need 

Q(ab) > 0, Q{ed) > 0, [Q{eb + ad)]^ < 49(a5)9(cd). 

Note that 

4%ab)%ed) - [^ieb + ad)]'^ = [^{cb - ad)]^ - \bc - ad\^. 

Minimizing the expression in (5.5) over 3?(w), we see that we also need 3?(c6 — ad) > if 
S(c(i) = 0, and 

(5.6) Q{ab) - [^(^^+^^)]^ + ^(^cb - ad)t + Q{cd)t^ > 0, t>0 

if ^{cd) > 0. Since the discriminant of this quadratic is 

[^{eb - ad)]^ + [^{eb + ad)]^ - 4^{ab)^{ed) = \be-ad\^>0, 
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if '^{cd) > 0, (5.6) is always true if be = ad, while if be ^ ad, (5.6) is equivalent to 3?(c5— ad) > 
0. Putting these observations together, noting that stability is not changed if the roles of b 
and c reversed, completes the proof. 

Proof of Theorem 5.1. We need to show that Qp{zi, Zn) ^ whenever ^{zi) > for 
i = 1, ...,n. To to so, fix zs, ...Zn with ^(zi) > for i = 3, ...,n, and write 

h{z, W) = Q{z, W, Z3, Zn). 

Then h is of the form considered in Proposition 5.2, and we must show that if (5.1) holds for 
a given a, b, e, d, then it also holds with b and c replaced by 

b{j>) = pb + {1 — p)e, e{p) — pe+ {1 — p)b. 

This follows from 

'^{ab{p)) — p'^{ab) + (1 — p)'^{ac), '^{ac{p)) — p'^{ac) + (1 — p)'^{ab), 
^{b{p)d) = pQ{bd) + (1 - p)^{ed), ^{e{p)d) = p^{ed) + (1 - p)^{bd), 

and 

3f?(6(p)c(p) - ad) = 3?(6c - ad) + p{l - p)\b - e\^ , 
b(p)e(p) — ad = (be — ad) + p(l — p){b — c)^, 

so that 

3?(6(p)c(p) — ad) — \b{p)e{p) — ad\ > "Ripe — ad) — \bc — ad\ 
by the triangle inequality. 

Corollary 5.3. Suppose the generating polynomial of the initial distribution of a symmetrie 
exelusion proeess on a finite set S is stable. Then the same is true at time t > 0. 

Proof. View the process in terms of stirrings. In other words, for each pair x,y E S interchange 
r]{x) and r]{y) at Poisson times at a certain rate. If stirrings are applied at only one pair of 
sites, the generating polynomial of the distribution at time t is of the form Qp given in the 
statement of Theorem 5, where Q is the generating polynomial of the initial distribution. 
For a general exclusion process, the distribution at time t can be obtained as a limit of that 
obtained by successively applying stirrings at different pairs of sites. 
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